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Abstract 

The magnetic moment /i of the p-meson is calculated in the framework of QCD sum 
rules in external fields. a<(-corrections are taken into account. Bare loop calculations 
(parton model) give: /ipart = 2.0. The nonperturbative effects (operators of dimension 6) 
result in negative contribution, whereas perturbative corrections - positive one. The final 
answer is = 1.8 ± 0.3. 
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1 Introduction 



QCD sum rule technique is an effective model independent tool in the study of var- 
ious characteristics of hadron states. This technique can be applied in particular to the 
investigation of the static properties of vector mesons. 

Namely, in the study of photoproduction, scattering, etc. of the vector mesons the 
vector dominance hypothesis is used. In the framework of this hypothesis it is supposed 
that the interaction of real or virtual photon with hadrons proceeds in such a way that 
the photon first transforms into vector mesons p, cj, 0, which then undergo interaction 
with hadrons. In the consistent lagrangian formulation of vector dominance model it 
is assumed (see j2] for review) that p-mesons are Yang-Mills vector bosons. In this 
hypothesis the p-meson magnetic moment is equal to 2 (in units eh/{2mpc)), at least if 
strong interaction is neglected. 

In paper [HI the p-meson formfactors were found at intermediate momentum transfer 
by QCD sum rules. By extrapolation of the p-meson magnetic formfactor to the point 

= (outside the applicability domain of the technique) it was found that the p-meson 
magnetic moment p is close to 2. However, this result can not be considered as final, and 
direct calculation of the magnetic moment is the topical problem. 

In our previous paper ^ the p-meson magnetic moment was calculated in the frame- 
work of QCD sum rules in external fields. This technique was developed in The 
operators of dimension 6 were taken into account. As for a^-corrections, they were ne- 
glected. However, since the discussed energy scale is approximately 1 GeV^, as-corrections 
may give appreciable contribution. 

This paper is devoted to the calculation of the cts-corrections to the p-meson magnetic 
moment. It is organized as follows. First of all, we are going to remind the reader the 
basic points of the sum rule construction (section 2). In the section 3 we discuss in detail 
the perturbative contribution to the sum rule. In the section 4 we present the nonper- 
turbative contribution. The results can be found in the section 5. The values obtained 
in some other approaches ( Dyson- Schwinger equation based models [7], jH], relativistic 
quantum mechanics jU], light cone QCD sum rules [10 J are also presented. 

In our paper u- and rf-quarks are considered as massless. 



2 Phenomenological part of the sum rule 

Let us consider the correlator of two vector currents in the external electromagnetic 

field: 



Here subscript F denotes the presence of the external electromagnetic field with strength 
Fpx and is the vector current with p-meson quantum numbers: jp = U'jpd. Its matrix 
element is 




(1) 
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In the limit of weak external field we consider only linear in Fpx terms in the correlator 
(1): 

n^, = nj, + iV4^n^.x-^x<x • (3) 

In order to find the p-meson magnetic moment we construct sum rule for the invariant 
function n(p^) at certain kinematical structure of ^^uxo^ (3)- 

The description of the kinematical structure choose and sum rule construction can be 
found in j3j. For the sake of consistency we present it here. 

The electromagnetic vertex of the p-meson has the following general form [3]: 

(p(p + g,e"')|jj|p(p,e")) = -e^'e^f ((2p + g)^^^,, - (p + g)^^^,, - + 



- 9x-ip)F2{-i^) + -^{p + q)pPcT{'2p + g)x^3(-g')) • (4) 

nip J 

In this expression j^' = e^^^Ai + edd'j^d is the electromagnetic current, , are u- 
and (i-quark charges and Fi , F2 , F3 are electric, magnetic and quadrupole formfactors 
correspondingly. The p-meson magnetic moment p is related to magnetic formfactor at 
g2 = 0: 

p = l + F2(0). (5) 
Using (2) and (4), we obtain for the fi\p) {p\jx\p) {p\ju\0)e^ transition: 



-^E(0|J.l/)(/lJ?l/)(p1J.|0)e, 

r,r' 
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= i^Y.^''p€^yv^x[{i'^P + l)x9po -{P + q)p9xa - Pa9px}Fii~q^) + (6) 

9p r,r' ^ 

+ (M<^ - 9x<T<lp)F2{-q^) + -^{p + q)pPa{'2p + g)xi^3(-g^)) • 

Here e-^— photon polarization, r, r' are the p-meson polarization indices. Let us consider 
in this expression linear in q„ terms. We sum over p-meson polarizations, retain the 
antisymmetric over x, a part, introduce the electromagnetic field F-^„ = i{e^q^ — e^q^ 
and obtain for (6): 

"^p _ 1 1 



2^"2^X-((^2 + -Fi)^{pu{Px9pa - Pa9^^x) - Pp{Px9ua " Pa9ux)) + 



11 \ 

+ 2^i^{pAPx9pa - Pa9px) + Pp{Px9va - Pa9ux)) + (^2 + Fi){g^^g^„ - gpa9vx)) ■ 
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Formfactor F3 does not give linear in g^- contribution. 

The p-meson magnetic moment /i is involved in the invariant functions at the kine- 

matical structures Pu{Px9^la - Pag^lx) ~ Pti{Px9yo - Pagvx) and g^x9ua - g^lagux■ We choose 
the structure 

PuiPxgfirr - Pag^x) - Pt^iPx9ua " Pa^i^x) ' ( 

In comparison with another one, (7) contains two additional powers of momentum in the 
numerator, which result in better convergence of the operator expansion series. 
Using equalities (5) and -Fi(O) = 1 one can see that ^2(0) + l/2Fi(0) = /i — 1/2. 
Thus, one should calculate the invariant function Il{p^) at the structure (7) in H^uxa- 
In order to obtain the phenomenological part of the sum rule, one should saturate 
the dispersion relation for Il{p^) by the contribution of physical states. We use the 
simplest model of physical spectrum, which contains the lowest resonance and continuum. 
Phenomenological representation of 11 has the form: 



n(p2) = J ds- 



Pis) 



p2y 
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Here dots mean the contributions of non-diagonal transitions (for example, 
(0|j^|p'^)(p*|j;^e^|p)(p|ji,|0), where p* is the excited state with the same quantum numbers 
as p), function f{s) represents continuum contribution and Sp is the continuum threshold 
for the p-meson. 

Retaining only the terms, which do not vanish after Borel transformation, we obtain: 



2g'j^ {m?p — p^Y m?p — p^ J [s — p^Y 
where C appears due to non-diagonal transitions. 



3 Perturbative contribution to the sum rule 

Now we need to calculate n(p2), basing on the operator product expansion in QCD. 
In this section we consider perturbative terms. 

The loop diagram (fig.l) corresponds to the operator of the lowest dimension 2. 



Figure 1: Loop diagram in the external electromagnetic field. Thick solid lines denote quarks, 
dashed line denotes the external field. 
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Figure 2: Perturbative corrections to the bare loop. Wave lines denote gluons. 



Its contribution to 11 (p^) is equal to 



ds 



2\2 ■ 



IGvr' J (s — p- 

The expressions for the quark propagator in the external electromagnetic field and in 
the external electromagnetic and soft gluon fields can be found in 

In order to improve the free-quark approximation, we calculate the as perturbation 
corrections corresponding to the diagrams in fig.2. We consider the integrals with respect 
to internal momenta in Z) = 4 — 2e dimensions. For the calculation of the diagrams in 
fig.2e and fig.2f the following formula from the paper [11^ is used: 



d^kd^q 



1 



(27r)2^ k^{p-kyq\q-ky{p-qy 



d^kd^q 



1 



(27r)2^ \k^{p - kYq^{k - qf k^{p - k)^q^{p ~ qf j ' 

In the limit e ^ these diagrams give the following contributions to the polarization 
operator: 

= ofi \ 9 I - + 22 - 87s ) (diagram in fig.2a) , 



a, / 4 
^^ = T^^^( 10 + 87ij 

2°'K-^P'^ V g 



(diagram in fig. 2b) 



a. 



( 14 + 87i5 



(diagrams in fig. 2c, d) 



«s /4 44 ^ 
A, = Af = :,^4^(- + — - 87^; 



2%V^e 3 



(diagrams in fig.2e,f) 



Here 7^; is the Euler constant. 
In the sum 

Aa + A, + + A, + + Af = 

the divergence 1/e cancels out. It is not by chance, since the diagrams in fig. 2 correspond 
to the physically measurable process, and their total contribution can not include diver- 
gent term. 

After adding the Os-correction, we obtain for the perturbative (and bare loop) contri- 
bution to the polarization operator: 

According to the quark-hadron duality, the continuum contribution in the interval 
of = —p^ from Sp to infinity is determined by term (9) in this interval. Therefore, 
function / in (8) is: / = -3/(167r2)(l + 10a,/(97r)). 



4 Contribution of dimension 6 operators 

The contribution of the operators of dimension 6 was found in j^]. Here we remind 
the corresponding formulas. 

The diagrams in fig. 1,2 correspond to the operator of the lowest dimension Fpx- The 
operators of dimension 4 are absent. As was shown in 0, operator q^Dpju — Dy'~^p)q 
has opposite with respect to the electromagnetic field C-parity and can not be induced 
by them, while operator e^up\(llblpD\q vanishes due to equation of motion for massless 
quarks. 




Figure 3: The diagrams, corresponding to the vacuum expectation value {qcrpxq)F{QQ)- Dots 
stand for permutations. 



Among the diagrams, corresponding to the vacuum expectation values of dimension 6 
operators (there are five such vacuum expectation values), the dominating contribution 
appears from the no loop diagrams with hard gluon exchange. In our case such diagrams 
contain the vacuum expectation value {qcTpxq) f{(1(i) (fig-3). Their contribution to the 
polarization operator is equal to 

9 p6 ^ ^ 

In this expression the quark condensate magnetic susceptibility x is negative. 

The contribution of the vacuum expectation value (G'"^G'"^)F^,y to H(p^) is determined 



6 



by the one-loop diagrams (fig. 4) and, therefore, has numerically small factor. Since these 
diagrams have infrared divergence, in order to evaluate them we introduce the cut-off over 
transversal momenta A and obtain: 

^ + -^"1 (U) 

36 ^TT ^V2AV QX^p^ p^^' ^ ^ 

As we shall see, (11) is really small in comparison with (10) for any reasonable A. 




Figure 4: The diagrams, corresponding to the vacuum expectation value 



The infrared divergence is probably canceled by the contribution of three other vacuum 
expectation values of dimension 6 operators: 

g{q{{G;,D^ - "d^G;,) - {G:,D, - 'D,G:,))^,eq)F , 

e,.px9{q{G;^D, + 'D,G;^)^a^eq)F , d^'' g' {G^.G^G^^ - GI,G\^G'^^,)) ^ . (12) 

Here is the covariant derivative, d^^'- are SU(3) structure constants, subscript F denotes 
the presence of the external field. Usually vacuum expectation values of such operators 
can be calculated by constructing the corresponding sum rules. But in our case this 
way is inapplicable because of their high dimension. However, the vacuum expectation 
values (12) are suppressed by the factor A*""^ {N^ is the color number) as compared with 
{q'^p\q)F{qq), and we disregard it. 

Collecting the expressions (9), (10), (11), we find the operator product expansion part 
of the sum rule: 

9, 3 lOttsx f ds 



^9 / 36^7r^^(2AV 6AV'^P^^' ^^^^ 



Now we have both parts of required sum rule. 



5 Results and discussion 

After Borel transformation 



^ ^ p2,„^oo n\ V dP^ ' 
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we equate the phenomenological (8) and operator product expansion (13) parts of the 
sum rule and obtain: 



Here C appears due to nondiagonal transitions. 

We use the following values of parameters: 
nip = 0.77 GeV^— the p-meson mass, 
gp/{4:n) = 1.27— the p-7 coupling constant, 
Sp = l.SGey^— the continuum threshold for p-meson, 
{{as/TT)G^) = 0.009 ± 0.007 Ge1/^-the gluon condensate [12], 
c,2(^g)2 ^ (0.28 ±0.09) X 10-2 Ge\/6 -the quark condensate 
X = —(5.7 ± 0.6) GeT^-^- the quark condensate magnetic susceptibility [13], 
AQCD = 0.2GeV. 

In (14) as(M2) = 47r/(9 ln(M2/AQ^25)). For the cut-off over transversal momenta A 
we take the value A^ = 0.8 GeV"^. 

As was mentioned in in order to obtain the p-meson magnetic moment in the 
parton model approximation one has to substitute the relation for Qp ^3] 



2 Sn'mi ■ " 



into the contribution of the bare loop and continuum (see (14)): 

_ 1 _ 3^. 

flpart ^ - ^^2^2 

It was obtained in this way: 

Impart 2 . 

This result agrees with the prediction of the vector dominance hypothesis. 

Now let us analyze the whole equation (14). In order to find the value of the magnetic 
moment, we approximate the right-hand side of (14) R{M'^) (see fig. 5) by a straight line 
in the interval 1.0 Ge^^ < JVP < 1.3 GeV^ and find it ordinate at zero Borel mass. 

Thus we obtain: 

p = 1.8. 

The operators of dimension 6 give less than 20% to this value. 

The contribution of the gluon condensate does not exceed 30% of the total contribution 
of dimension 6 operators, the contribution of the terms, which contain A^, is < 15% of the 
same quantity. That is why variation of A^ within the interval 0.6 GeV^ < A^ < 1.0 GeV"^ 
does not change the value of the magnetic moment. 

The variations of the values of the quark and gluon condensates within the given limits 
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Figure 5: The right-hand side of equation (14) i?(M^) as the function of IvP (sohd hue). The 
same without as-corrections (dashed hne). 

change the value of magnetic moment by < 10% each. 

The uncertainty in the value of the quark condensate magnetic susceptibility results 
in the error about few percent in the value of the magnetic moment. Variation of the 
continuum threshold for the p- meson Sp in the reasonable limits gives the same effect. 
Taking into account all uncertainties we obtain: 

/i = 1.8 ±0.3. (15) 

This is our final result. 

In jin] it was shown that the approximation procedure is correct (nonlinear terms can 
be safely neglected), when CM^/p <ti 1. In our case CM'^ / ^ ^ 0.2. 

Thus we see that in parton model (bare loop) approximation the p-meson magnetic 
moment fipart = 2. All accounted operator product expansion corrections are negative, 
i.e. nonperturbative interactions decrease this quantity. We obtain that the value of the 
magnetic moment (disregarding a^-corrections) is equal to 1.5 |3]. On the other hand, 
the a^-correction is quite appreciable (~ 0.3) in magnitude and positive (see fig.5). After 
taking into account both perturbative and nonperturbative terms, we obtain the value 
(15). It is consistent with the vector dominance model prediction again. 

The value of the p-meson magnetic moment was calculated in a number of papers 
within the Dyson- Schwinger equation based models. In ^ the value p = 2.69 was found. 
In jH] several results are compared, and the values of p lie between 2.5 and 3.0. Relativis- 
tic quantum mechanics model gives = 2.23 ± 0.13. 

Recently the value of the magnetic moment was found in the light cone QCD sum 
rule technique ^Oj- In this technique the expansion over the twist of the operators is 
performed, and all nonperturbative effects are encoded in the wave functions. In ^0] the 
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asymptotic form ip{u) = 1 of the twist 2 photon wave function ip{u) at the point u = 1/2 
was used, and it was obtained (disregarding Os-corrections): /i = 2.2 ± 0.2. This value 
and (15) do not contradict one another within the errors. 

The author is thankful to B.L. loffe for the statement of the problem and valuable 
discussions and to A.G. Oganesian for helpful discussions. 

The work is supported in part by grants INTAS 2000 Project 587 and RFFI 03-02- 
16209. 
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